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We study the dependence of the spectrum of four-dimensional SU(A0 gauge theories, where 
is the coefficient of the topological term in the Lagrangian, for N >3 and in the large-rV limit. 
We compute the 0(9 2 ) terms of the expansions around = of the string tension and the lowest 
glueball mass, respectively a(9) = a (l + S2& 2 + ...) and M(0) = M (l + g 2 9 2 + •••), where a 
and M are the values at 8 = 0. For this purpose we use numerical simulations of the Wilson 
lattice formulation of SU(N) gauge theories for N = 3,4,6. The 0(9 2 ) coefficients turn out to 
be very small for all N >3. For example, s 2 = -0.08(1) and g 2 = -0.06(2) for N = 3. Their 
absolute values decrease with increasing N. Our results are suggestive of a scenario in which 
the dependence in the string and glueball spectrum vanishes in the large-Af limit, at least for 
sufficiently small values of |0|. They support the general large-Af scaling arguments that indicate 
= 9 /N as the relevant Lagrangian parameter in the large-Af expansion. 
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1. Introduction 

Four-dimensional SU(A0 gauge theories have a nontrivial dependence on the angle that 
appears in the Euclidean Lagrangian as 

2 

(<7(jc) = ^i£iivpaF^ v (x)Fp' a (x) is the topological charge density). Indeed, the most plausible 
explanation of how the solution of the so-called U(1)a problem can be compatible with the l/N 
expansion requires a nontrivial dependence of the ground-state energy density F(6), 

F(0) = -^lnj [dA]exp (- J d d x^e] , (1.2) 

to leading order in l/N. The large-Af ground-state energy is expected to behave as [[I], [|] 

AF(d) = F(d) - F(0) = ^ d 2 + O (l/N 2 ) (1.3) 

for sufficiently small 0, i.e. < n. This has been supported by Monte Carlo simulations [^]. 
Indeed, numerical results for N = 3,4,6 are consistent with a scaling behavior around 0=0: 

/ '(0) = o- 2 AF(6) = ^C6 2 (\+b 2 d 2 + ...), C = Gc + c 2 /Af 2 + ..., b 2 = b X2 /N 2 + ... (1.4) 

where a is the string tension at = 0. C is the ratio %/c 2 where % = j d x(q(x)q(0)} is the 
topological susceptibility at = 0. Its large-Af limit Coo is ^ |J] Coo ~ 0.022. Estimates of c 2 and 
b 2 ,2 are [|] c 2 « 0.06, and b 2 , 2 « -0.2 (& 2 ~ -0.02 for SU(3) ). Eq. (Q can be recast as 

f(d)=N 2 f(d = d/N), /(0) = ^C0 2 (l+£ 2 2 + ...), (1-5) 

where Z? 2 = b 2t2 + 0(l/N ) = 0(1). This is consistent with 1/A^ scaling arguments, indicating 
6 = 8 /N as the relevant Lagrangian parameter in the large-A^ limit of the ground-state energy. 

The dependence of the spectrum is particularly interesting in the large-A^ limit, which may 
also be addressed by other approaches, such as AdS/CFT correspondence, see e.g. Ref. The 
analysis of glueballs using AdS/CFT suggests that the only effect of in the leading large-A^ limit 
is that the lowest glueball state becomes a mixed ++ /0 h state, but its mass does not change [g]. 

We present an exploratory study of dependence in the spectrum of SU(A0 gauge theories, 
using numerical simulations of the Wilson action. MC studies are very difficult in the presence of 
the complex valued term, defying simulation. We focus on relatively small 0, where one may 
expand observable values about 0=0. For the string tension and the lowest glueball mass we write 

a(0) = a(l+s 2 2 + ...), M(0)=M(l+£ 2 2 + ...) (1.6) 

where M is the ++ glueball mass at = 0. The dimensionless quantities s 2 and g 2 can be computed 
from correlators at = 0; they should approach constants as a — > 0, with 0(a 2 ) corrections. 

We present results for 4-d SU(/V) gauge theories with ,/V = 3,4,6. The O(0 2 ) coefficients turn 
out very small for all N >3, e.g., s 2 = —0.08(1) and g 2 = —0.06(2) for N = 3. Moreover, their 
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absolute values decrease with N. 0(6 2 ) terms are substantially smaller in dimensionless ratios 
such as Mjyfd and ratios of independent k strings, = o^jo. Our results suggest a large-A^ 
scenario in which 8 dependence in the string and glueball spectrum vanishes around 6=0. They 
are consistent with arguments indicating B = 6/N as the relevant parameter in the large-A^ limit. 
We also show a similar situation in the 2-d CP^ -1 models by an analysis of their l/N expansion. 

In Sec. | we outline the method to estimate 0(6 2 ) terms of the expansion in 6. The results of 
our numerical study are presented in Sec. [|. In Sec. |] we discuss the 6 dependence of 2-d CP^ -1 
models. Ref. [^] is a longer write-up of this work, with a more complete list of references. 



2. Numerical method 

2.1 Monte Carlo simulations 

In our simulations of lattice gauge theories, using the Wilson formulation, we employed the 
Cabibbo-Marinari algorithm to upgrade SU(/V) matrices by updating their SU(2) subgroups. This 
was done by alternating microcanonical over-relaxation and heat-bath steps, typically in a 4:1 ratio. 

The topological properties of fields defined on a lattice are strictly trivial. Physical topological 
properties are recovered in the continuum limit. Various techniques have been proposed to asso- 
ciate a topological charge Q to a lattice configuration; the most robust definition of Q uses the index 
of the overlap Dirac operator. However, due to the computational cost of fermionic methods and 
the need for very large statistics to measure correlations of Polyakov and plaquette operators with 
topological quantities, we used the simpler cooling method, implemented as in Ref. [J3]] . Compar- 
ison with a fermionic estimator shows good agreement for SU(3) |Q ^]. Moreover, the agreement 
among different methods is expected to improve for larger N. 

A severe form of critical slowing down affects the measurement of Q, posing a serious limi- 
tation for numerical studies, especially at large N. The available estimates of the autocorrelation 
time Zq for topological modes appear to increase as an exponential or a large power of the length 
scale pi ph. This dramatic effect has not been observed in plaquette-plaquette or Polyakov line 
correlations, suggesting an approximate decoupling between topological and nontopological ones, 
such as those determining confining properties and the glueball spectrum. But, as we shall see, 
such a decoupling is not complete. Therefore the strong critical slowing down that is observed in 
the topological sector will eventually affect also the measurements of nontopological quantities. 

2.2 The 0(6 2 ) coefficients of the 6 expansion 

Let us describe how to determine the 0(6 2 ) coefficients in Eq. (1.6). We first discuss the 



string tension; it can be determined from the torelon mass, i.e. the mass describing the large-time 
exponential decay of wall-wall correlations Gp of Polyakov lines. In the presence of a 6 term 

G P (t,9) = (A P (t))g = I[d j^}%!X^ . = I TrP (°;°) Tr^i,x 2 ;0, (2.1) 

P{x\,X2',t) is the Polyakov line of size L along the xj, direction. The time separation t is an integer 
multiple of the lattice spacing a: t = n t a. The correlation Gp can be expanded in 6 : 

G P (t, d) = Gf\t) + \d 2 Gf (/) + 0(d 4 ), (2.2) 
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where: d§\t) = (A P {t)) 6=Q , G [ p 2) (t) = -(A P (t)Q 2 ) e = + (Ap(t)) 6=Q (Q 2 )e=o ■ 
The correlation function Gp is expected to have a large-? exponential behavior 

G P (t,d)^B(d)e- E ^ t , (2.3) 

where E(8) is the energy of the lowest state, and B(8) is the overlap of the source with this state. 
If L is sufficiently large, the lowest-energy states should be those of a string-like spectrum. 

E(8) = a(8)L-n/(3L) (2.4) 



We expand the large-? behavior (23) of G(t, 8) as 

G P (t, 8) « B e- Eot [1 + 6 2 h(t) + ...] (2.5) 
where: B(d) = B Q + 8 2 B 2 + ..., E(8) = E + d 2 E 2 + h(t) = (B 2 /B ) - E 2 t . 



Comparing Eq. (2.5) with Eq. (2.2), we find that 



h(t)=G 2 2) (t)/(2G i0 Ht)) (2.6) 
Thus E 2 can be estimated from the difference: Ah(t) = hit) — hit +a). Indeed, lim Ah(t) = E 2 a. 

f— »oo 

Corrections are exponentially suppressed as exp[— (Eq — Eo)t], where Eq is the mass of the first 
excited state at 8 = 0. Assuming the free-string spectrum, Eq—Eq = An/L. Since we choose the 
lattice size L so that l a = JoL xi 3, (Eq* - E ) /E « 471//^ xs 1 .4. 



Finally, the dimensionless scaling coefficient s 2 of the 0(6 ) term in ( |1.6[ ) is obtained by 



E 2 

s 2 = ^ (2.7) 
oL 

s 2 is expected to approach a constant in the continuum limit, with 0(a 2 ) scaling corrections. 

An analogous procedure can be used for the lowest ++ glueball mass M(6). We employ 
wall-wall correlators of Wilson loops with up to 6 spatial links. We define: Ak(t) = k(t) —k(t + a), 
where k(t) is analogous to h(t), using glueball correlators. Then, g 2 in ( [IT^ ) is obtained by 

g 2 = lim Ak(t) (2.8) 
aM t-**> 

In order to improve the efficiency of the measurements we used smearing and blocking proce- 
dures to construct operators with better overlaps. Our implementation is described in Ref. [B]. 



3. Results 

Table [j] contains some information on our MC runs for Af = 3,4,6 on lattices L 3 x T. Since 
the coefficients of the 8 expansions are computed from connected correlation functions, and turn 
out to be quite small, high statistics is required to distinguish their estimates from zero: Our runs 
range from 9 to 25 million sweeps, with measures taken every 20-50 sweeps. This requirement 
represents a serious limitation to the possibility of performing runs for large lattices and in the 
continuum limit, especially for large N, due also to the severe critical slowing down. For all values 
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N 


j8 


lattice 


stat 


a 5 a 


aM ++ 


M ++/s/o 


3 


5.9 


12 3 x 18 


25M/20 


0.0664(6) 


0.80(1) 


3.09(4) 


3 


6.0 


16 3 x 36 


25M/40 


0.0470(3) 


0.70(1) 


3.23(4) 


4 


10.85 


12 3 x 18 


16M/50 


0.0646(6) 


0.76(1) 


2.99(5) 


6 


24.5 


8 3 x 12 


9M/50 


0.114(2) 


0.83(1) 


2.46(4) 



Table 1: Information on our MC simulations. Estimates of a are obtained using Eq. (2.4). 



of /3 considered, Tg satisfies Tq < 100 |^]. Furthermore, j8 values were chosen in the weak- 
coupling region (see Ref. [^] for a more detailed discussion of this point). The lattice size L was 
chosen so that l a = yfcL > 3 (see, e.g., Refs. [Q, ^]). Due to these limitations, in particular for 
N = 4,6, we could afford only one value of j8, so that no stringent checks of scaling could be 
performed. For this reason our study should be still considered as a first exploratory investigation. 
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Figure 1: Plot of Ah(t) and AJt(f) for N = 3 at J3 = 5.9 and /3 = 6.0. 

Figs. 1 and 2 show the results for the discrete differences Ah(t) and Ak(t) for N = 3 at j3 = 
5.9,6.0 and for N = 4,6 respectively. As expected, the signal degrades rapidly with increasing t. 
Anyway, they appear rather stable already for small values of t. In the case N = 3 the data at j3 = 6 
appear to approach the asymptotic behavior more rapidly than at j8 = 5.9. This should be due to 
the fact that more effective blocking can be applied when L = 16, rather than L = 12. 



We estimate S2 and g2 in (L6) from Ah(t) and Ak(t), taking the data at t/a = 2 in the N = 3,4 
runs, and at t/a = 1 for N = 6. In Table || we report the results. The estimates of S2 and g2 are 
small in all cases, and decrease with increasing N. For N = 3 the results at j8 = 5.9 and j8 = 6.0 
are consistent, supporting the expected scaling behavior. As final estimate one may consider 



s 2 = -0.08(1), g 2 = -0.06(2) for N = 3 



(3.1) 



One may also consider the the scaling ratio: M(G)/ \fc(d) = (M/y / a)(l +C2Q 2 + ■■■), where 
C2 = g2 — S2/I. Using the results of Table ||, we see that the 0(6 2 ) terms tend to cancel in the ratio. 
Indeed, we find c 2 = -0.02(2), -0.01(3), -0.01(2) respectively for N = 3,4,6. 

For N >3 there are additional independent ^-strings associated with representations of higher 
w-ality. One may consider the ratio Rk(6) = Ok(6)/o(6) = /?&(1 +^2 6 2 + •••)> where Ou is the 
^-string tension (see e.g. Refs. [Q, ^, [l(], |ll]]). For N = 4 there is one additional k string, 02 ; for 
N = 6 there are two. Our results for k > 1 strings are less stable. We obtained sufficiently precise 
results only for N = 4. They suggest a very weak 0-dependence in R2, i.e. \r%^\ < 0.02. 
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Figure 2: Plot of Ah(t) and Ak(t) for N = 4 at j3 = 10.85 and N = 6 at j5 = 24.5 . 



N 


j8 


S2 


82 


3 


5.9 


-0.077(8) 


-0.05(2) 


3 


6.0 


-0.077(15) 


-0.07(4) 


4 


10.85 


-0.057(10) 


-0.04(3) 


6 


24.5 


-0.025(5) 


0.006(15) 



Table 2: Results for S2 and g2, as derived from the discrete 
differences at t/a = 2 for Af = 3,4, and at t/a = 1 for = 6. 



In conclusion, the above results show 
that 0{Q 2 ) terms in the spectrum of SU(A0 
gauge theories are very small, especially 
in dimensionless ratios. Moreover, they 
decrease with increasing N, and the co- 
efficients do not appear to converge to 
a nonzero value. This suggests a sce- 



nario in which the 6 dependence of the spectrum disappears at large ,/V. General arguments indi- 
cate 6 = d/N as the relevant parameter, implying that O(0 2 ) coefficients in the spectrum should 
decrease as l/N 2 . This is roughly verified by our results, barring possible scaling corrections, es- 
pecially for N = 4, 6. For example, in the case of the string tension, S2 ~ S2.2/N 2 with S2.2 ~ —0.9 . 
Of course, further investigations are required to put this scenario on a firmer ground. 

Recent studies [ 12, [13] at finite temperature have shown that in the large-Af limit the topologi- 
cal properties remain substantially unchanged up to the first-order transition point. 



4. dependence in the two-dimensional CP N model 

Issues concerning the 6 dependence can also be discussed in two-dimensional CP^ -1 mod- 
els [|TJ], which present several features of QCD and, in addition, are amenable to a systematic 1 /N 
expansion around the large-Af saddle-point solution p4| , |15| ] . 

One may expand the ground state energy F(d) about 6 = 0. Defining a scaling quantity f(6): 

f(6) =M- 2 [F(d)-F(0)] = \cd 2 (\ + £ b 2n 6 2n ) (4.1) 

1 n=l 

M is the "zero momentum" mass at d = 0, and C is the ratio % /M 2 at d = 0, where x is the topolog- 
ical susceptibility. Within the 1 /N expansion: C = z/M 2 = \/(2nN) + 0{\/N 2 ). One obtains ^2n 
from correlation functions of q{x). The analysis of the 1/Af-expansion Feynman diagrams shows 
that b2n is suppressed as: Z?2« = 0{\/N 2n ). Thus the ground-state energy becomes 

f(e)=Nf(d = d/N), f(d) = ^Cd 2 (l + ^b 2n d 2n ), (4.2) 
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where C = NC and b 2n = N 2n b2 n are 0(1) in the large-V limit. Note the analogy with SU(AO gauge 
theories. The calculation of b 2n is rather cumbersome; we obtain 

b 2 = -27/5, b 4 = -1830/7. (4.3) 

Within the 1/7Y expansion one may also study the dependence of M on 6. We write: M(6) = 
M (l +m20 2 + ...). A diagrammatic analysis indicates that ni2 is suppressed as: ni2 = 0(1 /N 2 ). 
This confirms the arguments indicating B = 6/N as the relevant parameter in the large-TV limit. 
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